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ABSTRACT 

We evaluate the role of radiation pressure in the dynamics of H II regions. We first determine under 
what conditions radiation pressure is significant in comparison to gas pressure and show that, while 
radiation pressure is generally unimportant for H II regions driven by a handful of massive stars, it 
is dominant for the larger H n regions produced by the massive star clusters found near the Galactic 
center and in starburst environments. We then provide a solution for the problem of how H n regions 
expand when radiation pressure influences their behavior. Finally, we compare radiation-dominated 
H II regions to other sources of stellar feedback, and argue that H II regions are probably the primary 
mechanism for regulating the formation of massive star clusters. 

Subject headings: galaxies: star clusters — HII regions — ISM: bubbles — radiative transfer — stars: 
formation — stars: winds, outflows 



1. INTRODUCTION 

The dynamical expansion of a bubble of warm, ionized 
gas produced by a young star or star cluster is one of 
the classical problems of the interstellar medium. H II 
region dynamics are important in part because they play 
a do minant role in regulat i ng th e for mation of st a r clus - 
ters. IWilliams fe McKed (Il997f ) and ICarpenterl (|2000D 
argue based on observations that no more than ~ 10% of 
the mass in a giant molecular cloud can ever be incorpo- 
rated into stars, and Matzner] (|2002l ) and lKrumholz et alJ 
(|2006|) present GMC models that explain this inefficiency 
quantitatively in terms of the evaporation of GMC mass 
by H II regions. Nearby star clusters provide direct ev- 
idence for this phenomenon. For example in the Orion 
Nebula Cluster, the radiation of 9 1 Ori C launches an ion- 
ized wind from the molecular clo ud surface t hat carries 
a mass flux of - 1(T 2 M yr^ 1 (|O'delll200lJ ). sufficient 
to ablate a mass comparable to the ONC's stellar mass 
in ~ 10 5 yr. Any understanding of star formation rates 
and efficiencies, and of molecular cloud lifetimes, must 
therefore be based on an understanding of H II regions. 

Most work on this problem to date has focused on H n 
regions where conditions arc similar to those found in the 
Galaxy within several kpc of the Sun. The star clusters 
in this region are typically below 10 4 Mp, in mass (e.g. 
IWilliams fc McKedTl997t lLada fe 1. adl] |2( )(>■",:. and since 
a fully-sampled IMF at zero age prod uces S ~ 10 46 5 
ioniz ing photons s _1 per M Q of stars (|Krumholz et al.l 
l2006t) . they have luminosities of at most S ~ 10 50 ioniz- 
ing photons s _1 . For ionizing luminosities in this range, 
both numerical treatments and analytic estimates show 
that radiation pressure is generally sm all compared to 
gas pressure inside H II regions (e.g. [M athews 119691 : 
iGail fc SedlmavJl97ftlArthur et alJ200llHennevll2007M . 
Thus while radiation on dust grains may play an impor- 
tant rokj_rjroducmg srnafl holes in H n region centers 
(e.g. iMathewsl 119671 : Ihiouel 12002). the standard assump- 
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tion has been that radiation pressure is not significant 
for determining the dynamics the H n region as a whole. 
Treatments of H II regions under circumstances where ra- 
diation pressure is significant have generally been limited 
to h ydrostatic models that do not in clude any dynamics 
(e.g. iDopita et aLll2002t 120031 120061) or numerical mod- 
els of particular regions (e.g. iPellegrini 
lHarper-Clark fe Murravl[2009[ ) . The former provide little 
information on how radiation pressure affects gas motion, 
while the latter do not easily lead to general conclusions 
about when radiation pressure is important and how H n 
region dynamics are altered when it is. 

A re-assessment of the role of radiation pressure in H n 
region dynamics is timely because observations of star 
cluster formation are beginning to probe new regimes 
where earlier arguments minimizing the role of radiation 
pressure no longer apply. The importance of radiation 
pressure rises as the ionizing luminosity does, and in 



contrast to the values of at most S 



10 



50 



found 



in local H II regions, the Arches cluster near the Galactic 
center has an ion izing luminosity of S ~ 4 x 10 52 s _1 
(jFiger et al.ll2002l ). and some extragalactic clusters have 
even larger values. 

Moreover, regions with large ionizing luminosities also 
tend to have large virial and escape velocities, and in 
this case the conventional description of H n region ex- 
pansion driven by ionized g as pressure also breaks down. 
iMcCradv fc Graham! (p007h measure velocity dispersions 
of 10 — 30 km s _1 for the super star clusters in M82, 
and this probably represents a lower limit on the velocity 
dispersions and escape velocities in the clusters' parent 
molecular clouds. Since the sound speed in ionized gas 
is only 10 km s" 1 , an H II region driven by gas pressure 
could not have expanded within or driven mass out of the 
the clouds from which these star clusters formed. We are 
forced to conclude that either these clusters must have 
formed with an efficiency of nearly 100%, or that some 
mechanism other than ionized gas pressure is responsible 
for removing mass and limiting the efficiency. 
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Fig. 1. — Schematic representation of the blister (left) and 
embedded (right) H II region cases. The gray represents in the 
ambient cloud, and the black dot is the driving source. In the 
blister case the H II region is hemispherical, while in the embedded 
case it is spherical. 

While we cannot directly rule out the hypothesis of 
100% star formation efficiency for the M82 clusters, we 
can rule it out for similar large clusters in other galax- 
ies. For example, the young star clusters observed in 
the Antennae galaxies hay e com parable es c ape v eloci- 
ties (IWhitmore et al.l[1999l) . but iFall et all (|2005f) and 
IWhitmore et all (|2007l ) find that the star clusters in the 
Antennae show the same "infant mortality" phenomenon 
as those in the Milky Way: roughly 90% of all star clus- 
ters dissolve within ~ 10 Myr of formation, almost cer- 
tainly because the formation process operates with a low 
efficiency and removal of a majority of gas leaves the re- 
maining stars unbound. Clearly some mechanism must 
remove gas from these clusters as they form, and we ar- 
gue below that radiation-driven H II regions are the most 
natural explanation. 

Our approach to the problem is as follows. In §[2] we de- 
rive a condition for when radiation pressure is important, 
and we then give a solution to the idealized problem of an 
H II region expanding into an ambient medium including 
radiation pressure effects. In § [3] we discuss the contri- 
bution of trapped radiation to H n region dynamics. In 
§[4] we discuss the relative importance of radiation-driven 
H II regions and supernovae. Finally, we summarize in 

§E1 

2. THE EFFECTS OF RADIATION PRESSURE 
2.1. When is Radiation Pressure Significant? 

Consider a source of bolomctric luminosity L that pro- 
duces i onizing p h otons at a rate S, located at r — 0. Fol- 
lowing IMItznei] (|2002f ). we investigate two cases, which 
can be treated in parallel. The first is a region of neutral 
gas of density p = po{r/r Q )~ kp . The second is a region in 
which the density is p = po{r/ro)~ kp for x > and for 
x < 0. The former corresponds to the case of an "embed- 
ded" H II region that is completely surrounded by dense 
gas, and the latter to a "blister" H n region in which 
the driving source is at the edge of a dense cloud, and 
the ionized gas can escape freely We illustrate these two 
possible configurations in Figure [TJ We define e = 13-6 
cV as the threshold photon energy required to ionize a 
neutral hydrogen atom, and for convenience we define 
ip = L/(Seo) to be the ratio of the star's bolometric 
power to its ionizing power, counting only an energy eo 
per ionizing photon. For massive stars and clusters whose 
luminosity comes mostly from massive stars, ip ~ 1. 

We wish to determine how the gas moves in response to 
the radiation flux, and to understand under what circum- 
stances radiation plays an important role in determining 
gas motions. Following the usual procedure, we approx- 



imate that the ionized gas is isothermal at temperature 
Tji, and has a sound speed c\\ which is much larger than 
the sound speed in the neutral gas. If radiatio n pressure 
is negligible, we have the usual ISpitzerl fl978) solution: 
the ionized material expands due to its thermal pressure. 
This expansion sweeps the neutral gas into a thin shell, 
which contains most of the mass that was originally in- 
side the radius ni of the H II region. If radiation pres- 
sure is the dominant force acting on the gas then a fluid 
element at a distance r from the source undergoes a ra- 
diative acceleration a ra( j = J n v (y)Le~ Tv ^ j (knr 1 c) dv , 
where k„ (r) is the opacity of the fluid at radius r to pho- 
tons of frequency is, and T„(r) = J Q K v (r')p(r') dr' is the 
optical depth from the source to that point. Photons 
below the Lyman limit carry roughly half the radiative 
momentum, and since these are absorbed primarily by 
dust grains, K v (r)e~ T "( r ' jr 2 and thus the radiative ac- 
celeration is a decreasing function of r. The other half of 
the momentum is carried by photons above the Lyman 
limit, which can be absorbed by either H I or dust grains 
- see Appendix lAl If dust absorption dominates, then ra- 
diative acceleration falls with radius as for lower-energy 
photons. If H I absorption dominates, the acceleration is 
proportional to the recombination rate, which is either 
flat (if the gas density is uniform) or again declines with 
radius (if radiative acceleration causes gas to pile up near 
the shell edge). Thus, the total radiative acceleration is 
always largest closest to the source, and again material 
will be swept into a thin shell of radius rri- The interior 
of this shell will be optically thin. 

Thus, after a rapid initial expansion phase, the dynam- 
ics of the gas reduce to the problem of computing the 
dynamics of the th i n she ll that bounds the H II region. 
Following [Matznerl (|2002h . we can solve this problem by 
writing down the momentum equation for the shell: 



dt 



(M sh r n ) = A sh <{ pn[cjj + mii(wh - ni)] + 



/trap L 



(1) 

where M s h and Ajh are the shell's mass and area, and 
pu and mii are the density and velocity of the gas im- 
mediately interior to it. The shell area and mass are 
A sh = (4, 2)7rrf T and M sh = (4, 2)irrf{p/3, where p{r) = 
[3/(3 — k p )]p (r/r )~ k p is the mean density inside the 
spherical or hemispherical region of radius r in the initial 
cloud, and the values in parentheses refer to the cases for 
a (spherical, hemispherical) H n region. The first term 
on the right-hand side represents gas pressure, while the 
second represents radiation pressure. Note that in writ- 
ing this equation we have implicitly assumed that all the 
radiation force is applied at the thin shell, rather than in 
the H II region interior. This is certainly a good approx- 
imation when radiation pressure is dominant, since, as 
discussed above, the interior of the shell will be cleared 
by radiation pressure, and thus all photons will be ab- 
sorbed in or near the shell. When gas pressure dominates 
and the shell interior is of uniform density, the rate of 
momentum deposition by ionizing photons matches the 
recombination rate. Since this is uniform, so the mean 
radius at which momentum is deposited is 3/4 of the shell 
radius. Moreover, non-ionizing photons, which carry half 
the total momentum, still deposit all their momentum in 
the shell. Thus our approximation that all the momen- 
tum is deposited in the shell is a good one. 
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The quantity /trap in equation {T]) represents the fac- 
tor by which the radiation pressure force is enhanced 
by trapping of energy within the expanding shell. If 
/trap = 0, then the shell is optically thin and all stel- 
lar photons escape without depositing any momentum. 
Given the high color temperature of the emitting stars 
and the large opacity of the gas clouds where massive 
clusters form, this is not realistic. A value /trap = 1 
corresponds to every photon emitted by the stars being 
absorbed once in the shell and depositing its momentum 
there, then promptly escaping. If there is more than 
one interaction per photon then /trap could potentially 
be much larger than unity. In this case the velocity to 
which the shell accelerates will be limited by the rate at 
which the stars supply energy rather than momentum. 

Trapping can happen in three ways. First, some 
fraction of the stellar radiation will go i nto accelerat- 
ing lin e driven winds off stellar surfaces (|Castor et al.l 
Il975bh . and the expanding wind will collide with the 
slower-moving shell. This will produce some transfer 
of momentum, which could be larg e if the shocked gas 
becomes trapped inside the shell ([Castor et al.l Il975at 
IWeaver et al.lll977f ). Second, if the shell is sufficiently op- 
tically thick to long-wavelength radiation, then ultravio- 
let and visible photons that are absorbed by dust grains 
in the shell and re-radiated at infrared wavelengths may 
remain trapped in the shell and interact more times be- 
fore finally escaping. Third, Lyman a photons that are 
produced by recombinations in the shell or in the H II 
region interior may undergo many resonant interactions 
before escaping. We defer a discussion of these trapping 
mechanisms until § |3l and for now we simply assert the 
result from that section: /t ra p is always likely to be of or- 
der a few. For this reason, we choose to leave it as a free 
parameter of constant value, for which we adopt a fidu- 
cial value /trap = 2 when we wish to perform numerical 
evaluations. 

We can characterize when radiation pressure is sig- 
nificant by examining the limiting cases of gas- and 
radiation-pressure dominated flows. In the g as-pre ssure 
dominated case, we have the usual ISpitzerl <| 19781 ) H n 
region solution. Once expansion of the H II region be- 
comes subsonic with respect to the ionized gas, the H n 
region interior approaches a uniform density, and ioniza- 
tion balance requires that 



4 3 ( PiiCii V 



(2) 



where ru, Xn, and Cn are the radius, temperature, and 
sound speed of the ionized region, qb is the case-B recom- 
bination coefficient, is a dimcnsionless number that ac- 
counts for absorption of ionizing photons by dust grains 
and for free electrons provided by elements other than 
hydrogen, and we have adopted the usual on-the-spot ap- 
proximation. If He is singly-ionized and 27% of photons 
are absorbed by dust rather than g expected for gas 
pressure-d ominated H H regions with Milky Way dust-to- 
gas ratios dMcKee fc Williamsl[T997T) . then <j> = 0.73. We 
discuss the value of 0, and of dust absorption generally, 
in more detail in Appendix [AJ Note that equation ^ 
holds approximately even in the c ase of a blister- type 
hemispherical H n region. Following Matzner (2002|), we 
consider in that in the case of an embedded, spherical 
H II region that uu <C cn, while for a blister- type one 



uu w 2cn- Thus the gas pressure term on the right-hand 
side of equation |T]) becomes 



pn[cn + uii(un - hi)} ~ (1, 2) 



3S(j) ksTn 
4ttq:r r 3 / 2 



(3) 



In the limiting case of a radiation-pressure dominated 
flow, the radiation pressure term is simply L/(47rr 2 T c). 
Since the radiation- and gas-pressure terms have differ- 
ent radial dependences, we can calculate a characteristic 
radius for which they are equal: 



12(l,4)7r0 \k B Tu 
(9.2,2.3) x lCT 2 ^ pc, 



J trap 9 



(4) 
(5) 



where the numerical evaluation is for our fiducial param- 
eters Tn = 7000 K, = 0.73, / trap = 2, tp = 1, and 
a B = 3.46 x 10~ 13 cm 3 s" 1 , and 5 49 = S/10 49 s" 1 . Since 
radiation forces vary with radius as r^ 2 , radiation dom- 
inates at smaller radii and gas pressure at larger radii. 

It is useful to compare this to the Stromgren radius at 
which gas pressure-driven expansion begins in the case of 
negligible radiation pressure force. Setting pu in equa- 
tion j2]) equal to the mean density p( r st,o) inside rg tj o, 
we find 



p{rst,o)4i 



\47raB / 

and computing the ratio of r cri to this we find 



c 



»*st,o 



trap 



1,2)/cb?iic/ \36tt4> 2 iiumn 



(6.2,1.6) x 10- 2 (n H ,2S 4 9) 2/3 , 



(8) 



where pu = 0.61 is the mean molecular weight in the 
fully ionized gas (so cn = 9.74 km s _1 ) and 7Th.2 = 
p(^St,o)/(100/zmH) is the mean density of H nuclei in 
units of 10 2 cm -3 , and p = 1.4 is the atomic mass per H 
nucleus for gas of standard cosmic composition. 

Thus we see that for single OB stars, S49 ~ 1, expand- 
ing into Galactic molecular clouds, Sh,2 ~ 1, radiation 
force is negligible once the H n region reaches a tenth 
of a pc in size (r cri < 0.1 pc), and is often negligible as 
soon as the ionized region has finished its initial rapid 
expansion to the Stromgren radius (£ < 1). Thereafter 
the usual gas pressure-driven expansion solution applies. 
However, we reach a very different conclusion if we con- 
sider the formation of very massive clusters in dense en- 
vironments. In Table [T] we list properties for a sample of 
massive star clusters in the Milky Way, M82, the Anten- 
nae, and NGC 5253. We plot S49 versus 7Ih,2 for these 
objects in Figure [5J As the plot shows, these clusters 
have C in the range ~ 1 — 10 4 , indicating that they go 
from the border between radiation- and gas-dominated 
to completely radiation-dominated. The characteristic 
radii where gas pressure becomes comparable to radia- 
tion pressure run from ~ 1 — 100 pc, which is generally 
larger than the physical size of the cluster in question. 
The exception is the Orion Nebular Cluster (the inverted 
triangle), by far the smallest cluster shown in Figure [2l 
This has r crl = 0.06 pc (assuming a blister-type H 11 re- 
gion, which is observed), considerably smaller than the 
0.8 pc radius of the cluster. 
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References. — 1 = IHillenbrand fc Hartmannl (Tl998l1 . 2 = I Williams fc McKed (IT997I1. 3 = 
IFiger et dl999l'>. 4 = IKim et al.l 112000 1 , 5 = ITurner etaLI (120001 , 6 = IMcCradv fc Graham! 
I|2007l ). 7 = IGilbert fc Graham! H2007l i 

Note. — Col. (2): Stellar mass. Col. (3): Half-light radius, except for Antennae objects, 
where we use 1/2 of FWHM. Col. (4): Density, computed as n H ,2 = [3M/(47i\R 3 )]/(100/wn// ). 
Col. (5): Ionizing luminosity. Col. (6-8): r ch , and r sta n, computed with fiducial parameters 
in the blister case with k p = (equations |51 [71 and !2H . 

a M82 objects are super star clusters, whose identifiers correspond to those given in 
IMcCradv fc Grahaml ((2CI07h b An ionizing luminosity for this object has not been reported in 
the literature, so we estimate it by t aking S49 = 10 — 2,5 (M/Mq), the value for a fully-sampled 
zero-age IMF (Krumholz ct al. 2006). c Antennae objects are emission line clusters, whose iden- 
tifiers correspond to those given in Gilbert & Graham (2007) d This radius is almost certainly 
an overestimate of the radius at formation, since the A ntennae emission line clusters have un- 
dergone significant expansion l|Gilbert fc G raham 2007]). Correcting for this effect would raise 
^H,2j which in turn would increase f and decrease r sta n. 



We should add a caution here, which is that for the 
clusters described in Table Q] and plotted in Figure [2 we 
have computed r c h, C: an d fstaii using the full ionizing 
luminosity of the cluster. This is not correct very early 
in the expansion process, since the initial Stromgren ra- 
dius around each massive star is so small that it may en- 
close at most a few of its neighbors. Since this will lower 
the ionizing luminosity compared to our value, the very 
early expansion could be gas-driven. However, as the ex- 
panding shells overlap, the incorporation of more stars 
into their interiors will rapidly convert the expansion to 

a radiatively-driven one; since ( cx S% , the expanding 
shell in a cluster for which we have computed C, = 100 
reaches C > 1 an d becomes radiation-dominated when 



the shell includes only 4.6% of the cluster luminosity. 
Since luminosity tends to be strongly centrally concen- 
trated, this will occur well before the shell includes this 
fraction of the cloud volume. Thus radiation takes over 
very early in the expansion process for those clusters for 
which we have computed values of £ 3> 1 . 



2.2. The Dynamics of Expansion with Radiation 
Pressure 

To study the dynamics of an expanding H 11 region 
with significant radiation pressure, it is convenient to 
non-dimcnsionalize the equation of motion ((T|). We let 
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Fig. 2. — 549 versus ?Th,2 for a sample of massive star clusters. 
The value of r cll shown on the right axis is derived from equation 
JSJ, and the lines of constant f (dotted) and r sta n (dashed) are 
derived using equations J7]l and 12H . All are for the blister case. 
Values of r ata ii indicated in the figure are in pc. The data shown 
are described in Table [TJ emission line clusters in the Antennae 
Galaxies (red circles), super star clusters (SSCs) in M82 (blue tri- 
angles), the SSC in NGC 5253 (purple star), and the Arches (green 
square) and Orion Nebula (brown inverted triangle) Clusters in the 
Milky Way. 

x = r/r c h and r = t/t c h., where 



tch 



1 47T p r Q p c ^- kp 

3 kp ftra,pL 



(9) 



With this change of variables equation (fTJ) becomes 



d 



3-fcp d 
x u -xn )=1 



„l/2 



(10) 



The first term on the right hand side represents radiation 
pressure, and the second represents gas pressure. This 
equation is not exactly correct near in = 1, because in 
writing the gas pressure term we have implicitly assumed 
that the density inside the H 11 region is uniform and 
has the value given by ionization balance, equation @. 

This assumption cannot be precisely true when x\\ 1, 
because while radiation pressure is significant it will ex- 
ert a force on material in the H 11 region interior that 
will make the density higher toward the shell wall than 
near the H 11 region center, while equation assumes 
uniform density. Nonetheless, equation (|10[) represents 
a reasonable approximation that becomes exact every- 
where except near x\\ = 1. 
At early times, when xu <C 1, we can drop the gas prcs- 
1/2 

sure term x :l on the right hand side, and the resulting 
equation admits the similarity solution 



^Il.rad 



4 - k B 



l/(4-fc„) 



(11) 



The dynamics of this solution can be understood by 
noting that since the momentum of the shell equals 
the radiant momentum modified by the trapping factor 
(M sh f 11 = ft ra pLt/c), the shell's kinetic energy M sh f u /2 
is a very small fraction, / t 2 ap fn/(2c), of the total ra- 
diated energy Lt. In fact, if /trap ~ 1, then the en- 
ergy of the shell's motion approximately matches the 
energy of the photons currently crossing its interior: 
M sh f I 2 I /2 a (r)/2)Lr n /c, where 77 = r„t/r„ = 2/(4- k p ). 



The inefficiency of direct radiation driving is related to 
the low value of /trap, an issue we return to in § [3] 

If we were to drop the radiation pressure term, we 
would have the usual similarity solution for gas pressure 
expansion 



3*1,1 



(7 - 2fc p ) 2 
4(9 - 2k p ) 



2/(7-2fc p ) 



(12) 



Although the exact solution will approach this value 
when xu 3> 1, the existence of a characteristic scale r c h 
implies that there is no true similarity solution that in- 
cludes both the radiation- and gas-driven phases. When 
radiation pressure is significant it provides an extra boost 
of momentum in the phase before gas pressure takes over, 
accelerating the expansion relative to the standard sim- 
ilarity solution. This breaks the self-similarity of the 
gas-driven expansion phase, although the expansion ap- 
proaches self-similarity as r — ► 00 and the extra momen- 
tum provided by the radiation becomes small compared 
to that input by the gas. In the absence of a similarity 
solution, however, it is trivial to integrate equation (|10p 
numerically. We do so subject to the boundary condi- 
tion that xu and dxn/dr approach the values that cor- 
respond to the similarity solution ([IT]) as t — > 0, and plot 
the result in Figure [3] for some sample values of k p . For 
comparison, we also show the pure radiation and pure 
gas similarity solutions. The true solution may be rea- 
sonably approximated by an appropriately weighted sum 
between the two. The expression 



3-II,approx 
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(7-fe p )/2\ 2/(7 k ^ 
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(13) 



is accurate to better than 5% for k p = — 1. 



2.3. H II Region Stalling 

As discussed in § [TJ one of the reasons that radiation- 
driven expansion is of interest is that gas-driven expan- 
sion will fail in regions where the ambient velocity dis- 
persion and escape velocity exceed the ionized gas sound 
speed (e.g. TlVIatznerl [20021 iKrumholz et al.ll2006f) . The 
expansion velocity is fu = v^dxii/dr), where 



'cli 



v ch = — = v/(378)C (2fep - 3)/4 cn. 

tch 



(14) 



The expansion rate therefore drops to the ionized gas 
sound speed when dxu/dr = C, lyZ ^ 2klD ^ A / >/(3, 8). Since 
we have solved for dxu/dr numerically, it is trivial to 
numerically invert this equation for a given £ to obtain 
the radius and time at which the expansion becomes sub- 
sonic. However, it is more illuminating to consider the 
radiation-dominated case ^ ^> 1 , because in that case 
the similarity solution pip applies, and the resulting an- 
alytic form for dxu/dr allows us to solve for the radius 
at which the subsonic transition occurs in closed form. 
In this limit, expansion becomes subsonic at 



*£sub 
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2(3,8) 
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To give some idea of this in physical units, for k p = the 
radius at which the expansion becomes subsonic is 



r sub = (0.90,0.90)(^4 9 /n H ,2) 1/2 



pc. 



(16) 
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Fig. 3. — Dimcnsionlcss radius xn (top row) and velocity dxn/dr (bottom row) versus time r. In the left column we show the numerical 
solutions for k p = 0, 0.5, and 1. In the middle we show a comparison of the num eric al solution and with the similarity solutions for pure 
radiation (da shed line, equation 1110 . pure gas pressure (dot- dashed line, equation 1120. and an appr oximation that combines them (circles, 
equation 1 1 3 D . In the right column we show the error in our approximate solution (equation 1136 , defined as the difference between the 
approximate and numerical solutions, normalized to the numerical solution. 



Clearly the transition to subsonic expansion will not hap- 
pen until an H II region has swept up a significant frac- 
tion of the gas in the protocluster, and it is therefore 
likely that radiation-driven H n regions can expel gas 
from clouds where gas-driven ones could not. 

Alternately, we can ask when the H n region expansion 
velocity becomes comparable to the velocity dispersion 
within the parent cloud. This should be a reasonable es- 
timate of when the expansion will stall. To the extent 
that the parent cloud is in virial balance, this will also a 
good estimate of at what radius the expansion velocity 
will be reduced to the point where gravity can prevent 
expansion. 1 We write the velocity dispersion a of mate- 
rial within a radius r of the cloud center as 



a vir GM{< r) 



5r 



1/2 
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a v i r Gp(rst,o>c h C "x 



1/2 



(17) 



(18) 



where M(< r) is the mass within a radius r of the ori- 
gin, x = r/r c h, and a v ir = 5a 2 r/\G M (< r)] is the virial 
ratio of the material inside radius r ( Bertold i fc McKeel 
[19921) . The factor (4,2) appears because the gas is spher- 
ical in the embedded case and hemi-sphcrical in the blis- 
ter case. The velocity dispersion can be either larger or 
smaller than en depending on the density and radius: for 
k p = 0, we have er/cn = 0.011(a v ir«H,2) 1 ^ 2 (?'/pc). The 

1 One might expect that comparing radiation forces to gravita- 
tional ones would not produce a characteristic radius, since both 
vary as 1/r 2 . However, this argument is only valid if the grav- 
itational force is dominated by a central point mass. Forming 
star clusters, however, are dominated by gas mass, and for a gas 
mass distribution p oc r~ k p, the gravitational force varies as r 1 ~ k p. 
Thus it is possible to obtain a characteristic radius at which gravity 
balances radiation. 



expansion velocity of the shell v c h(dxn/ dr) is equal to 
this velocity at the dimcnsionlcss radius xu for which 



dxii 



(4, 2)ira vh Gp{r S t,o)r; 
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-4fe p )/4 (2- 
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(19) 

As with the problem of determining when the expansion 
becomes subsonic, it is simple to find a numerical solution 
to this equation for a given choice of nu,2 and S49, but it 
is more illuminating to consider the radiation-dominated 
case, for which we can find a solution analytically. If 
we set xu equal to the value given by equation pip for 
the radiation-dominated phase, then the dimensionless 
radius x s t a ii that satisfies equation (| 1Q[) is 



X s ta\l 



(45, 180)4 



2tt(4 - k p )a vil Gp(r S t,oK h 



(4fe p -3)/2 



1 l/(4-2fe„) 

(20) 



Putting this in dimensional terms, for k p — and a, 
1 the stalling radius is 

r s t a ii = (8.9, 10.6)r2H2 /25, 49 4 pc. 



(21) 



Numerical evaluations of r s taii for our sample clusters are 
given in Tablc[TJ For most of the clusters shown, the stall 
radius is comparable to or larger than the total cluster ra- 
dius. This demonstrates that, unlike gas pressure-driven 
H 11 regions, radiation pressure-driven ones can poten- 
tially expel gas from the large clouds that form mas- 
sive clusters. Wc show lines of constant r sta n (computed 
for the radiation-dominated case), with observed clusters 
over-plotted, in Figured 

3. RADIATIVE TRAPPING 

Our derivation of the shell motion in § [2] assumes that 
during the radiation-dominated phase shells are momen- 
tum driven, in the sense that the force exerted on the 
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shell is the radiation momentum flux multiplied by at 
most a modest enhancement factor /t rap - However, if 
stellar radiation energy becomes trapped in the bubble 
interior, then the full amount of bolomctric energy emit- 
ted by the stars can be used to drive expansion, and 
/trap 3> 1- Here we estimate /trap considering three pos- 
sible sources of trapping: 



/trap 1 ^ /trap,w /trap,IR ~l~ /trap,Lyo 



(22) 



where the 1 represents absorption of the direct radia- 
tion and the remaining three factors on the right are the 
contributions to trapping due to stellar winds, infrared 
photons, and Lya photons. 

3.1. Stellar Winds 

Massive stars emit strong winds, during their main se- 
quence evolution and afterward. The observational and 
theoretical status of these winds is somewhat uncertain, 
due to the effects of line blanketing (which affects the 
mappi ng between spectral type and effective temper- 
ature; iMartins et all l2002t 120051 : iRepolust et ail l2004f ) 
and of wind dumpiness (which affects the inferred mo- 
mentum flux) . Nevert heless the wind-luminosity relation 
(|Kudritzki et al.ll!999h. which appears to re main valid af- 
ter these corrections (jRepolust et al.ll2004f ) . implies that 
a young cluster massive enough to sample the entire ini- 
tial mass function emits about one half as much force 
in winds as in starlight, M w v w = f w L/c with /„, w 0.5. 
The ratio f w reflects the flux-averaged optical depth of 
line overlap in the winds from hot stars, and so indicates 
a transfer of momentum from photons to wind. 2 

In a dense cluster, the winds emitted by different stars 
rapidly collide with each other and shock-heat up to tem- 
peratures of ~ 10 7 K. The hot, shocked gas will then 
accele rate to supersonic speeds in a thermally-driven 
wind (|Cant6 et a.1.1 l2000h . creating a bubble in the sur- 
roun ding ~ 10 4 K gas. Models of stellar wi nd bub- 
bles dCastor et al.lll975al IKoo fe McKelll992allbf l typi- 
cally assume that motions are efficiently energy-driven, 
rather than momentum-driven, so long as the wind is fast 
enough for its stopping shock to be non-radiative. An 
energy-driven bubble expands much more violently than 
the momentum-driven solutions presented above: the ki- 
netic energy of a radiation-driven shell is fu/(f w v w ) ~ 
10~ 2 times smaller than the wind energy. Thus our solu- 
tion is va lid only if the expansion is mo mentum-driven. 
However, lHarper-Clark fc Murray! (|2009f ) point out that 
wind energy is easily lost if the shell is porous so the 
shocked wind gas can escape. Indeed, such escape is man- 
dated by observations showing that the x-ray luminosity 
of large bubbles is far below what would be expected if 
all of the stellar wind energy remained trapped. 

We use a modified version of lHarper-Clark fc Murrayl s 
analysis to quantitatively estimate what this implies for 
/trap,w Suppose that the shell of material swept up the 
by the expanding H II region covers a fraction Cf of 
the sky as seen from the driving cluster; for blister case 
clearly Cf < 1/2, while for the embedded case Cf could 
take any value between and 1 depending on the shell's 

2 In optically thick Wolf-Rayet winds and LBV outbursts, by 
contrast, the wind momentum can exceed L/c; however the stellar 
energy budget is smaller after core H burning, and we assume the 
momentum budget is as well. 



porosity. Stars in the cluster launch winds with a mass 
flux M w at a velocity v w , and the fast wind material 
shock- heats to a temperature T x Tjj. The hot gas 
exerts a pressure Px = Pxc x on the shell, where p\ is 
the density of the shocked material within the shell and 
cx is its isothermal sound speed, so 



/fa 



pxc 2 x 



(23) 



/ trap,w — Tit* 2 \ ' 

The pressure of hot gas is set by the balance between 
the injection and escape of mass and thermal energy. 
Both are injected directly by stellar winds and escape 
through holes in the shell, but additional mass is incor- 
porated through the interaction of hot gas with the dense 
shell. We assume this stripping occurs through hydrody- 
namical ablation rather than thermal evaporation, partly 
because magnetic fields inhibit thermal conduction. Sim- 
ilarly, stellar photons only heat gas to temperature Tn, 
not to Tx, so they do not contribute; heating of gas 
to Tx by x-ray photons is relatively insignificant. The 
equations for mass and energy conservation of the hot 
gas are 



M X =M W - 47rrf T (l - Cf)p x c x + M, 



abl 



(24) 



E x =L w -4mri l {l-Cf) -p x 



- inr^CfPxru- (25) 
where middle terms in both equations assume the hot 
fluid escapes through holes at speed cx and has its aver- 
age properties as it does so. Leakage ensures that mass 
and energy do not accumulate, implying that the LHS 
of both equations is effectively zero; furthermore the last 
term in E w , which represents adiabatic work, is negligible 
in the presence of holes because Cffu *C (1 — Cf)c x . 

The term M a bi is the rate at which ionized gas at tem- 
perature Tn ablates off the inner edge of the shell and 
heats up to T x through its interaction with the wind. 
The ablation process is a complex one, and the rate of 
ablation depends on the thermal behavior of the shocked 
wind. When the cooling time within the mixing layer 
is long, as it tends to be when the hot gas originates 
in a fast stellar wind, mass a blation occurs i n the "jet- 
limited" regime described bv lCanto fc Ragal (|1991f ). In 
this case the ablation rate is limited by the rate at which 
the wind can supply momentum to a thin mixing layer 
between the two fluids, and these authors estimate that 
the ablation rate per unit area is 0.09/9ttc 2 t / (2cx), where 
pu is the density of the ionized layer abutting the mixing 
layer. Assuming this process occurs over the entire shell, 
the total ablation rate is 

Since the ionized gas is in pressure balance with the hot 
gas, puCjj = pxc x - Inserting this and L w = M w v^/2 
into equations (j2~4"|) and (f2"5"|) and omitting the negligible 
terms, 

M w = Anr 2 n pxc x [(1 - C f ) - 0.045C/] (27) 

(28) 



2 ^207rr 2 llP xCx(l~Cf). 
The solution is 



M w v w - 



pxc x = 



M u ,v w /(4 7 rr I 2 I ) 



[5(1-C/)(1 - 1.045C/)] 
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(29) 
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so that, using (1 - C/)(l - 1.045C/) ~ (1 - 1.02C/) 2 



firs 



fu 



a/5(1 - 1-02C/)' 



(30) 



The divergence as Cf approaches unity is not real, as our 
neglect of adiabatic losses and of accumulation within 
the shell are incorrect when the holes close up. Likewise 
values of /trap.w less than f w are not realistic, because 
the wind force is always present; our error in this limit is 
to assume that the wind energy is thermalized, when in 
fact it remains mostly kinetic if the shell is mostly holes. 

Since f w ~ 0.5, this implies that /t ra p,w — 0.22/(1 — 
Cf) ~ 1 unless we are considering the case of an em- 
bedded H II region with an extremely non-porous shell, 
1 — Cf <C 1. Such low porosity is implausible given the 
turbulent, clumpy nature of the interstellar medium and 
the fact that pressure-driven shocks tend to run down 
density gradients and "blow out" . Even if one started 
with a perfectly uniform ISM, the expanding shell is sub- 
ject to IVishniacl (|l983f ) instability of a pressure-driven 
slab; moreover if the wind caused the expanding shell 
to accelerate then the Rayleigh- Taylor instability would 
spontaneously create holes in the shell, reducing Cf be- 
low unity. Thus we conclude that /trap.w is at most a 
few, and is likely to be small than order unity. Our 
co nclusion is consisten t with the numerical simulations 
of iTenorio-Tagle et al.l (|2006f ) . who also find that, in a 
non-uniform medium, the bulk of the mass around a 
young star cluster is not swept up into the thermal wind 
driven by ~ 10 7 K gas. Instead, that gas escapes rapidly 
through the porous shell, while the bulk of the mass ex- 
pands more slowly (see their Figure 9). 

The above calculation depends somewhat on our es- 
timate of the term M a bi, which is uncertain for several 
reasons: because the physics of ablation is not well un- 
derstood, because the ablating area could be very differ- 
ent from 47rC/rj I , for instance if the shell's structure is 
more interesting than a broken sphere, and because other 
mechanisms like thermal evaporation, photoevaporation, 
and cloud disruption can all inject mass. Similarly, the 
density of the ablating gas could be higher than indicated 
by our pressure-balance argument, since photoevapo- 
ratio n might compress the gas (e.g. iBertoldi &: McKeel 
|1990() . Given these uncertainties, one might consider cx 
better constrained than M a bi - for instance, an upper 
limit on the X-ray luminosity implies a lower limit on 
cx- Taking cx as given and solving equations (|24|) and 
for Px, we find 



JV. 



trap.w 



5(1 - C f ) c x 



(31) 



so trapping is important when (1 — Cf)cx/v w < f w /5 ~ 
0.1. 

Protostcllar winds represent a separate, brief, but very 
intense phase which deserve separate mention. Being 
magnetically launched, they are strongly collimated and 
may far exceed the photon momentum; moreover the en- 
tire stellar population generates them. Rather than as- 
sess th em directly we appeal to the treatment bv lMatznerl 
(|2007| ). who found that protostellar winds are very signif- 
icant in the formation of the Pleiades and Orion Nebula 
clusters, but quite insignificant in the formation of the 
Arches or more massive clusters. 



3.2. Trapped Infrared Photons 

If the expanding shell traps the infrared light emitted 
within it, then its momentum will exceed that of the 
driving starlight by a factor 

Pm 



,/*trap,IR 



L/i^cY 



(32) 



where Ptr is the pressure of the trapped IR radiation 
field inside the shell. We first consider the highly ideal- 
ized case of a uniform, non-porous shell, Cf = 1, which 
provides an upper limit on /t ra p,iR; before treating leak- 
age through holes. 

When a uniform shell is optically thick to its own 
thermal emission, it radiates from a photosphere above 
which the flux-averaged optical depth is 2/3. The emis- 
sion is characterized by the shell's effective tempera- 
ture, which satisfies 47rrjl l (7sB7 1 c 4 ff sh = L, and the flux- 
averaged mass opacity above the photosphere is approx- 
imately the Planck mean ftp(T ff. sri ), implying a column 
S p h ~ 2/[3Kp(T c ff iS h)] above the photosphere. Shells 
with E sri = M s i 1 /(47r7'j I ) < £ pri are optically thin to re- 
processed light. Thick shells, with £ sri 3> £ p h, can be 
treated in the diffusion approximation (dPm/ kr(T) = 
— o'sbTsS w ith -Pir = aT(£) 4 /3 and kr the Rosse- 
land mean mass opacity), below the photosphere. The 
solution for the pressure within an optically thick shell is 
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(33) 



where P(Pr) = J^ m dP{ R /K R (T'), and is the in- 

verse function of J r . Evalua ting equation (|33|) for the 
IWeingartner fc Drainel (|2001f ) standard dust model "A" 
for Ry = 5.5, we find that 



fu ■ap,IR : 



yr 3 



132 

Toff, si 



V-l-92 



72 



P 



cll'.sh 



1.71 



-2/3 



(34) 

to good accuracy, where £ s hcii is understood to be in 
gcm~ 2 . This result shows that radiation trapping can 
be quite significant, for instance when £ s heii ^ 1 and 
Toff,sh > 60 K, as is expected around a luminous young 
cluster; however this estimate is unrealistically high when 
radiation can leak away. 

In the more realistic case where leakage is important, 
we can compute Pir in a manner analogous to our calcu- 
lation of Px, by balancing the rate at which energy the 
stars inside the shell add energy against the rate at which 
it leaks out. Our treatment h ere is a simplified vers ion of 
that given in Appendix D of iMcKee fc Tanl (|2008[ ). We 
limit our attention to the case where the shell is optically 
thick (Sghoii ^> £ p h(Xoff iS h)) 011 average: in this limit we 
can compute the energy density of the trapped radiation 
field when the shell is porous, Cf < 1, simply by treating 
the shell as a perfectly opaque sphere with holes in it. In 
this case energy balance requires that 

L = (1 - Cj)L + 3^(1 - C f )P m . (35) 

Here the left hand side is the rate at which the stars in 
the cluster inject radiant energy. The first term on the 
right hand side is the rate at which this beamed radia- 
tion field escapes the shell without interacting, while the 
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second term represents the rate at which the reprocessed 
infrared radiation field leaks out. In writing this term, 
we have approximated that the IR radiation is isotropic, 
so the radiation flux escaping through a hole in the shell 
is given by the expression for the flux through an in- 
finitesimal pinhole in an oven: F = 'ZcPnx/A. If the holes 
have finite size the true flux will be higher because the 
radiation field in the vicinity of the hole will be beamed 
outward. Solving this equation for Pm gives 



Pm = - 



C 



f 



L 



3 \1 -Cfy /./jjl 

and substituting this into equation (|32[) gives 



,A;rap,IR 



1 - c 



(36) 



(37) 



As with /trap,w, this is an upper limit. The true value 
will be smaller because radiation leaks out of the shell 
by diffusion at a finite rate even when C/ = 1, because 
we have assumed that the internal radiation field has 
had an infinite amount of time to reach its steady state 
value, and because we have neglected the enhanced rate 
of energy loss through holes of finite size compared to 
those of infinitesimal size. 

Even with this overestimate, our analysis shows find 
that for realistic values of Cf < 1/2, the ratio of the 
force provided by trapped IR to that provided by the di- 
rect radiation field is /trap,m 1- We further note that, 
even if we started with a perfectly uniform ISM and did 
not have any winds to punch holes in the shell, numer- 
ical simulations show that a trapped IR radiation field 
is in itself sufficient to induce a Rayleig h- Taylor-like in- 
stabi lity that punches holes in the shell (|Krumholz et al.l 
[2009h . This will ensure that Cf is always well below 1, 
so that /trap.m is no more than a few. 

3.3. Trapped Lyman a Photons 

Even blown-out sections of the shell that arc transpar- 
ent to infrared radiation may still be opaque to Lyman a 
photons, since they interact via a resonant process that 
produces a far larger cross section than absorption of IR 
photons by dust grains. As a result, our analysis based 
on Cf does not apply to Ly a radiation . However, as 
pointed out by iHennev fc Arthur! ([1998) , the influence 
of Ly a trapping is severely limited by the presence of 
dust grains. Absorption by dust grains in the H n region 
interior or its bounding shell converts Ly a radiation to 
IR, which can then escape through the porous shell as 
discussed in the previous section. This process limits the 
energy density and pres sure that can build up in Lyman 
a. We refer readers to IHennev fc Arthur! for a detailed 
treatment of this problem, and here simply quote the re- 
sult: the Lyman a radiation pressure saturates at a value 
Phya ~ 0.06/ciiiCjjcrjj), where er d! o = <7d/5 X 10~ 22 cm" 2 
and (id is the dust cross section per H nucleus for Ly a 
photons. Milky Way gas has (Td,o ~ 1- We therefore have 
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(38) 



Thus during the radiation-pressure dominated phase, 
xn < 1, trapped Ly a radiation is negligible in compar- 
ison to direct radiation pressure unless the dust is less 
opaque than typical Milky Way gas by a factor of > 20. 



Once gas pressure begins to dominate, x\\ > 1, trapped 
Ly a radiation may exert more pressure than direct stel- 
lar radiation, but in this case it will still be negligible in 
comparison to gas pressure. For this reason we can sim- 
ply set /trap,Lya ~ without making a significant error 
in any phase of the expansion. 

4. RADIATION-DRIVEN H o REGIONS VERSUS 
SUPERNOVAE 

We have argued that radiation pressure dominates the 
dynamics of H n regions in massive clusters, but this by 
itself does not establish that radiation-dominated H n re- 
gions are the dominant source of feedback in star clusters. 
The primary alternative model is that supernovae are 
the main form of feedb ack (e.g. iParmentier et al .1 [20081 : 
iBaumgardt et all [2008). However, we can show on the 
basis of time scales that supernovae can play only a lim- 
ited role. The crossing time in an object (either star 
cluster or gas clo ud) of mass M and surface density £ is 
(|Tan et al.ll2006f) 



tr. 



0.95 / M 



1/1 



-1/2^1/4^3/4 Myr ^ 



(39) 



wher e a v ; r is the object's virial ratio (Bertold i fc McKed 
fl992l) . Mi = A//(10 4 M Q ), and S = E/(l g cm" 2 ). 
Both observations and theory indicate that gas clumps 
that are forming massive star clusters have T,q ^ 1 
dMcKee fc Tad 12001 iKrumholz fc McKed l200l . and a 
virial ratio S> 1 would preclude the possibility of form- 
ing a gravitationally bound cluster. Thus we expect 
t CI ~ 10 5 yr unless the cluster mass exceeds 10 6 Mq. 
This theoretical argument is in good agreement with ob- 
served star cluster crossing times in cases where those can 
be me asured directly. For example. iMcCradv fc Graham! 
(|200l find crossing times of ~ 10 5 yr for 21 clusters in 
M82. Thus we conclude that the crossing time in pro- 
toclusters is at least an order of magnitude shorter than 
the ?s 3 Myr required for the first supernova explosion. 
Numerical studies of star cluster formation show that, in 
the absence of a source of energy to drive the turbulence 
(such as stellar feedback), star formation generally pro- 
ceeds to ef ficiencies of order un ity within a few crossing 
times (e.g. iBonnell et al.ll2003[ ). Thus supernovae could 
be an effective form of feedback only if some other, differ- 
ent energy source, most likely a different form of stellar 
feedback, inhibited star formation and kept the proto- 
cluster from collapsing for ~ 10 crossing times. While 
there is some observatio nal evidence for such an extended 
star formation history (jTan et al.l 120061 : iHuff fc StahTerl 
I20061 ). even then supernovae would at most remove the 
last vestiges of gas from the cluster after most of the 
gas had been consumed by star formation or expelled by 
other mechanisms. 

This argument argument does not by itself rule out the 
possibility that supernovae from neighboring star clusters 
that formed earlier might act as a source of delayed feed- 
back. However, this idea also has timcscalc problems. 
In Orion, for example, there have been o nly 10 — 20 SN 
explosions in the last 12 Myr (|Ballvll2008D . a rate of 1 - 2 
per Myr. Thus the typical cluster-forming cloud is not 
hit by a SN blast wave within its 0.1 Myr crossing time. 
Moreover, even if a blast wave did hit a cluster- forming 
cloud, any shock strong enough to affect a cloud signifi- 
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cantly tends to destroy it completely, ra ther than leaving 
it int act but exciting turbulent motions (|Nakamura et al.l 
[20M) . Thus if SN from neighboring clusters were a dom- 
inant form of feedback, then most cluster-forming gas 
cloud would their lives being shredded by a blast wave. 
This is inconsistent with the observed morphologies of 
disrupting clouds, which typically show no evidence for 
interaction with a strong external shock. 

5. SUMMARY 

We examine the role of radiation pressure in deter- 
mining the dynamics of expanding H 11 regions around 
young star clusters. We derive a characteristic number 
(equation [7]) that describes the importance of radiation 
pressure in H 11 region expansion, and we show that for 
very massive star clusters such as the Arches or the large 
clusters found in starburst galaxies, this number is al- 
ways 3> 1, indicating that radiation pressure dominates 
at least the initial phases of expansion. For such clus- 
ters, we derive a solution describing the expansion of a 
shell of material driven by a combination of radiation and 
gas pressure. Our solution shows that radiation pressure 
is capable of driving expansion at speeds considerably 
larger than the ionized gas sound speed, and therefore 
radiation-driven H 11 regions can be an effective feed- 
back mechanism in clusters where a high escape speed 



renders gas-driven H 11 regions ineffectual. 

Radiation-driven H 11 regions are likely to be the domi- 
nant feedback mechanism that regulates the formation of 
massive star clusters. Unlike supcrnovae, they begin to 
operate immediately upon formation of the stars, rather 
than with a delay of many dynamical times. While stel- 
lar wind bubbles or trapped infrared radiation fields are 
potentially more effective than direct stellar photons as a 
source of feedback, we argue that they are unlikely to be 
dominant because they are sapped by leakage of energy 
out of the expanding shell. Future models of feedback in 
massive star cluster formation must therefore be careful 
to take radiation pressure-driven H 11 regions into ac- 
count. 
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APPENDIX 

A. THE ROLE OF DUST ABSORPTION 

For most of the time that an H 11 region is expanding, dust within it will be able to absorb photons. Assuming 
that all of a source's radiant energy is eventually absorbed by dust grains (even if ionizing photons are first absorbed 
by hydrogen atoms and are eventually down-converted to photons below the Lyman limit), the equilibrium dust 
temperature at a distance r from the source is given implicitly by 

4^2 =/ Q A BAT)dv, (Al) 

where Q A is dust the absorption efficiency, which depends on grain properties. Dust can surv ive wherever T is below the 
sublimation temperature T s ; for the most refractory grains IWolfire fc CassinellH ((1986, 1987) estimate j Q A B V {T S ) dv = 
1.24 x 10 7 erg cm -2 s _1 , so at least some dust exists at all radii greater than the sublimation radius 

r s = 3.8 x lO-V 72 -^ 2 pc (A2) 

Even for the most luminous sources we consider, S*4g ~ 10 4 , this is much less than the cluster radius, r c h, or r sta ii- 
However, dust absorption does not affect the dynamics of the radiation pressure-dominated phase directly, since the 
dust and gas are likely to be well-coupled by magnetic fields. The momentum transferred by radiation to the gas 
therefore does not depend on whether photons are absorbed by dust grains or hydrogen atoms. However, dust grains 
do play a role when gas pressure becomes important, because they modify the ionized gas density pu and thus change 
the gas pressure; this dependence is captured through our parameter 4>. A full discussion of th e behavior of dust in side 
H II regions is beyond the scope of this work, and we refer readers to lDopita et all (|2002h and lArthur et all (|2004h . In 
t his appendix we mak e some simple estimates to show that our use of a constant <f> is unlik ely to lead to serious error. 

lArthur et alj (|2004f ) show that the value of </> (which is equivalent to the parameter y 3 in lArthur et all s treatment) 
is determined by the dust optical depth T d through the ionized gas; for typical gas pressure-dominated H 11 regions 
whos e metallicity is not vastly greater than Solar, tj, < 1 and <fi > 0.5 is appropriate, consistent with our value of 0.73 
from lMcKee fc Williams! (fl997T ). To see how this might vary for the radiation-dominated case, it is useful to estimate 
t c i in terms of the ionization parameter 

U = — , (A3) 
cn 

where F is the ionizing photon flux into the recombining gas and n is its density. Ionization-recombination balance 
requires <fiF = a.Bn 2 h, where h is the thickness of the recombining gas layer and n is its number density, so the dust 
optical depth is 

r d = nha d = ^-U ~ 9O0W, (A4) 
a B 
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where ad is the dust cross section per H nucleus to ionizing photons. The numerical evaluation uses ad — 1.0 x 10 
cm 2 , appropriate for Solar mctallicity. 

Since <j> < 1, this result shows that Td < 1 as long as hi < 0.01. When gas pressure dominates, n is constant 
and an increase in ionizing luminosity increases hi. However, once radiation pressure becomes significant, there is 
a countervailing effect: radiation pressure piles up gas against the inner wall of t he expanding shell, increasing its 
density and thus driving hi down. Under the assumption of hydrostatic equilibrium, iDopita et all (|2002h find that this 
effect ensures that hi never rises above ~ 0.01. While our expanding shell is of course not in hydrostatic equilibrium, 
as long as the expansion is slow compared to the signal speed within the shell, quasi-hydrostatic equilibrium should 
be established in the shell interior, and thus we also expect a maximum value of hi ~ 0.01. This ensures that dust 
absorption never dominates over ionized gas absorption, and we can safely adopt a constant, order-unity value of <\>. 

However, we note that even if this argument is incorrect, we only make a significant error in following the dynamics 
when xii ~ 1, i-e. when gas and radiation pressure are nearly of equal importance. That is because our choice <fi = 0.73 
is appropriate when gas pressure dominates (xn 3> 1), and does not affect the dynamics significantly when radiation 
pressure dominates (xn 1). Since our treatment of the crossover regime is very approximate in any event, inaccuracy 
in estimating <p in this regime would not further degrade our accuracy. 
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